Angle-resolved photoemission spectra in the cuprates from the rf-density wave theory 
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Angle-resolved photoemission spectra present two challenges for the d-density wave (DDW) theory 
of the pseudogap state of the cuprates: (1) hole pockets near (it/2, tt/2) are not observed, in apparent 
contradiction with the assumption of translational symmetry breaking, and (2) there are no well- 
defined quasiparticles at the antinodal points, in contradiction with the predictions of mean-field 
theory of this broken symmetry state. Here, we show how these puzzles can be resolved. 
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At first glance, the ci-density wave (DDW) proposal for 
the pseudogap state of the cupratesi seems to naturally 
explain the principal anomaly in photoemission spectra 
in this state: the existence of a gap with d x i_ y i symmetry 
without superconductivity. However, since DDW order 
breaks translational symmetry, thereby splitting the Bril- 
louin zone into two magnetic Brillouin zones, the Fermi 
surface in the first magnetic Brillouin zone should be du- 
plicated in the second magnetic Brillouin zone. Thus, 
the Fermi surface consists of hole pockets, which is of 
importance in understanding a number of experiments, 
such as superfluid density 2 , Hall numbed, etc. However, 
in angle-resolved photoemission spectroscopy (ARPES) 
in hole doped cuprates, Fermi arcs - not hole pockets - 
are observed*. There is spectral weight in the first, but 
not the second magnetic zone. In this paper, we show 
from a careful analysis that Fermi arcs rather than hole 
pockets are indeed the consequences of the DDW theory 
in ARPES. 

A second important aspect of the proposal of a broken 
symmetry state, even one of an unusual variety, is that 
it is expected to support electronic quasiparticles which 
are essentially Fermi-liquid like, as they are in a BCS 
superconductor. However, from ARPES in underdoped 
samples no peak is observed at the antinodal points in 
the normal state, but one appears in the superconducting 
state upon cooling. 4 This observation also finds a natu- 
ral explanation within our theory^ We show that the 
antinodal quasiparticles, being relatively high-energy ex- 
citations, decay by creating particle-hole pairs along the 
Fermi arcs in the DDW state. In contrast, in the d-wave 
superconducting state (DSC), or in the coexisting DDW 
and DSC state, the Fermi arcs shrink to points, and the 
decay rate is considerably reduced, resulting in a peak 
in the spectral function. This reduction is bolstered by 
the suppression of the decay matrix element by the su- 
perconducting coherence factors. 

The explanation discussed here involves interaction be- 
tween quasiparticles, whose absolute magnitude is set 
by a reasonable Hubbard-like interaction of magnitude 
1.5 eV, but the precise magnitude is of not much conse- 
quence. There may be other sources of broadening of the 
quasiparticle peak, including fluctuation effects, bilayer 
splitting, fractionalization, etc., which we do not address 
here. We merely wish to point out that within the sim- 



plest mean field picture of DDW, there are no obvious 
puzzles. 

In our mean field analysis, and indeed in many theo- 
ries, the nodal quasipartcles, or excitations at the Fermi 
arcs, should in principle be sharp, which is not entirely in 
keeping with ARPES, although a fairly well defined peak 
is observed both above and below the superconducting 
transition temperature, T c . It remains to be seen if the 
present experimental situation changes with time or not. 

To establish our notation, we begin with a brief sum- 
mary of the mean field theory of DDW. The Hamilto- 
nian TL can be simply written in the first magnetic zone 
by introducing the Pauli matrices o~ x , o~ z , the identity 
matrix I, a row vector ^ a = (c^ Q , —ic^+Q a ), and its 
Hermitian adjoint. The electron destruction operators of 
momentum k and spin a are Ck, Q and the momentum 
Q = (-7T, 7r). Thus, K, = TL — fjM is given by 
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W k a x ] * k>a , (1) 



Here [L is the chemical potential and M is the number 
operator. Note that K, is complex Hermitian, reflecting 
broken time reversal symmetry. In the first magnetic 
zone, it is convenient to define ej^ = \ [et ± cic+q] , where 
£k is the electronic band structure. A standard Bogoli- 
ubov transformation diagonalizes the Hamiltonian, but 
since I commutes with a x and cr z , (e£ — fx) can not enter 
the coherence factors, which are 
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where E k = \f(e^] 
trade places as k 



The coherence factors must 
k + Q, which is a consistency check 

li) . The 

energy eigenvalues are: 



as to why they cannot be functions of (ej 



E±=e+±E(k). 
The DDW gap is assumed to take the form: 
W (T) 



(cos k x — cos k y ) 



(3) 



(4) 



The electron spectral function in a crystal need not be 
invariant under translation by a reciprocal lattice vector. 
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In fact, it is weighted by the Fourier transform of the 
relevant Wannier orbitals. If the Wannier orbitals are 6- 
functions, the spectral weight is the same in all Brillouin 
zones. On the other hand, if the Wannier orbital is spread 
out spatially, then the spectral weight in higher Brillouin 
zones will be very small, and I(u>, k+G) <C I(e, k), where 
7(e,k) is the angle-resolved photoemission intensity. In 
the DDW state, the unit cell has been doubled. The 
coherence factors itk, ^k tell us how the two sites within 
the unit cell are superposed, so «k/ u k plays the role of the 
Wannier function. The corresponding spectral function 
in the DDW state is: 

^l=ul5(u J -E++^+vl6(u-E k +v) (5) 

Consider \l < 0, the case of hole doping, such that the 
chemical potential lies entirely in the valence band, so 
that E+ - yL > 0. Then the ARPES intensity is 

I(u,k) ocn F (uj)vl5(uj-E^+fj) (6) 

Since i^k+cj = Uk, the photoemission intensity in the first 
and second magnetic zones differ only by these coherence 
factors: 

I(u,k) oc n F (iv)vlS (u; - E k + fi) , (7) 
7(cj,k + Q) oc n F (uj)ul5(uj-E k + fi) . (8) 

For k in the first magnetic zone (i.e. for k + Q in the 
second magnetic zone), itk vanishes when Wk vanishes. 
In other words, the photoemission intensity in the sec- 
ond magnetic zone vanishes along the diagonals. For 
wavevectors close to the diagonals, the intensity goes as 
W£. Thus, the outer section of the hole pockets will have 
small or even vanishing spectral weight, and may not be 
detected in ARPES experiments. The spectral weight at 
a typical point on the outer part of a hole pocket will 
depend on various details, including the band structure, 
the precise angular dependence of the DDW gap, etc. 
A commonly-used model for the band structure is given 

by 

= At' cos k x cos k y , (9) 
= — 2i(cos k x + cos ky). (10) 

A generic parameter set is t = 0.3 eV, t'/t = 0.3, 
fi = —0.3 eV; with this set of parameters, the doping 
level is 14.3% The Fermi surface with a typical value of 
Wo(0) = 0.06 eV consists of four hole pockets as shown 
in Fig. ^ The corresponding i>£ appearing in the photoe- 
mission intensity is shown Fig. [3 It is clear that only half 
of the hole pockets will be visible in the ARPES spectra, 
resulting in Fermi arcs, despite the fact that the actual 
Fermi surface consists of hole pockets^ The results are 
similar for the second set of band structure parameters. 

We now turn to a discussion of the lifetime of a quasi- 
particle at the antinodal region k* close to (tt, Q), where 
the free electron Fermi surface crosses the band edge. 
The equation that determines k* is obtained by solving 



(0 ,K) (71 , 7C) 
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FIG. 1: The Fermi surface for Wo{0) = 0.06 eV. The band 
structure parameters are defined in the text. 




FIG. 2: The coherence factor u£. The parameters are the 
same as in Fig. Q 

for k y . Antinodal quasiparticles at k* have an energy 
very close Wo, the maximum of the DDW gap. Hence, 
they can scatter into a nearby wavevector while creating 
a particle-hole pair near the Fermi arcs (the inner section 
of the hole pockets). This is very different from the situa- 
tion in the ci-wave superconducting state, where there are 
only Fermi points, not arcs, as a result of which, there is 
very little phase space for low-energy particle-hole pairs. 
Secondly, the density of states is enhanced at the gap 
edge, resulting in an abundance of available phase space 
into which the quasiparticle can be scattered. In the su- 
perconducting state, this density of states enhancement 
is cancelled by coherence factors. These coherence fac- 
tors reflect the fact that the quasiparticles are neutral, 
so they are only weakly scattered by interactions which 
are coupled to charge. 

We will set up the lifetime calculation in full generality, 
assuming that both DDW and DSC order parameters are 
present, and then vary the size of the DSC order param- 
eter. In order to more easily compare with experimental 
results, we will assume mean-field-like temperature de- 
pendence for the DSC gap so that we can display our 
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results as a temperature-dependent decay rate. 

Consider an initial quasiparticle state of momen- 
tum ki in the antinodal region (to be precise, k* 
defined above) and of energy £ki> where £ k = 

yj (£L(k) — i±f + |A(k)| 2 . In the pseudogap state, 



where the d-wave superconducting order parameter A = 
0, £ k = S-(k) — [m. Suppose that this initial state decays 
into a final state of energy £k 2 +£k 3 + £k 4 • In lowest order 
perturbation theory, the decay rate for such a process is 



1 

n 



2tt 



'k 2 k 3 k 4 



|M klk2k3k4 | 2 (2 7 r) 2 (5(k 1 -k 2 -k3-k4) ( 5(fk 1 -fk 2 -^k3-fk 4 )[l-/(fk 2 )][l-/(fk3)][l-/^k 4 )], (12) 



r 



where Jk 2 k 3 k 4 =JllWWWWF> M kik 2 k 3 k 4 is a matrix 
element, and /(£ k )is the Fermi function. We have in 
mind a situation in which k 2 is close to ki and k3 , k4 
are close to the zone diagonal, but we will perform the 
integrals over the full Brillouin zone. 

There is a second contribution to the decay rate, I/V2, 
resulting from scattering off thermally-excited quasipar- 
ticles. The corresponding expression involves a different 
matrix element N klk2kski an d the quasiparticle at mo- 
mentum k 2 is thermally excited with probability /(£ k2 ). 
In all other respects the equation is the same as Eq. IT51 
except that the energy and momentum conserving 5- 
functions must be changed accordingly. The total decay 
rate is the sum 1/r = 1/Yi + l/r 2 . 

The matrix elements M klk2kaki and iV klk;2k3kl will de- 
pend on the form of the interaction between quasiparti- 
cles and also on the coherence factors. If we choose a 
Hubbard-like density-density interaction, Apj(q)/?|(— q), 
between the original electrons, the coherence factors are 
extremely complicated in the coexisting DDW and DSC 
state, and the multidimensional numerical integrations 
involved in calculating the scattering rates become next 
to impossible to carry out. To obtain upper bounds, we 
replace them by their maximum values. In the DDW 
state, we expect the coherence factors to be rather tame, 
but in the state with both DSC and DDW orders, they 
will suppress the decay rate as in a pure superconducting 
state. Thus, such an approximation will underestimate 
the difference between the decay rates in the pseudogap 
and the underdoped superconducting states. We call this 
interaction, treated with this approximation, model A. 

In order to capture the effect of the coherence factors 
in the mixed DDW and DSC state, we consider a model 
interaction. Since we are only concerned with the inter- 
action between the quasiparticles in the valence band, 
we choose the interaction to be (f2 is the volume of the 
system.) 



v = 75 E ^{^rVft^-qi. 



(13) 



q.k.k' 



where ip ka creates a quasiparticle in the valence band of 
the DDW state. We have ignored the temperature and 
momentum dependence of the interaction, because we 



are primarily interested in temperatures much lower than 
the DDW transition temperature, where the temperature 
dependence of the DDW gap should be weak. Moreover, 
we merely wish to demonstrate how the development of 
superconductivity affects the lifetime, so we also neglect 
the momentum space structure of the interaction. For 
this interaction, which we call model B, the coherence 
factors are equal to unity for the DDW order alone, but 
are non-trivial in the state with both orders as a result of 
the coherence factors associated with superconductivity 
in the mixed state. We can view the mixed state as DSC 
developing on top of DDW. Thus, the coherence factors 
for this interaction can be read off from the BCS theory 
and the matrix elements are 

M k k . k . k; = \[-V kl V k2 U k3 V^ + V kl U k2 V k3 Vu4 

- U kl V k . 2 U ks C/ k4 + U kl Vk 2 Vk 3 C/ k4 ] (14) 

and 

-^ki,k 2 ,k 3 ,k 4 =A[— ^ki Vk 2 Vk 3 Vk 4 — U kl U k2 U k3 U ki 

- C/ kl ^k 2 ^k 3 C/k 4 - V kl U k2 U k3 V ki 
+ U kl V k2 U k3 V ki - V kl U k2 U k3 V ki 

- U kl V k2 V k3 U ki + V kl U k2 V k3 U ki 

where 



(15) 



1 f E_(k)-// 

2 ^ £(k) 



(16) 



This form of the interaction allows us capture the differ- 
ence between the matrix elements in the pseudogap and 
superconducting states. 

Our results are displayed in Fig. |3 where the decay 
rates are plotted against temperature. In these plots we 
have kept the total gap fixed and equal to 0.06 eV, 
while assuming a mean-field temperature dependence for 
the superconducting gap 



\AT) = \ l w[l-^ r 



1/2 



(17) 



with A (0) = 0.03 eV and T c = 60 K. This implicitly 
defines the temperature dependence of the DDW gap, 
which is weak close to T c , as noted earlier. 
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FIG. 3: The lifetime of a quasipartcle in the antinodal region 
at wave vector k*, as defined in the text, plotted as a function 
of temperature. The open symblos are for without coherence 
factors and the solid symbols are with coherence factors. The 
square symbols correspond to 2 meV broadening of the energy 
conserving 5-function and the circles to 1 meV broadening. 
The electron-electron interaction parameter A = 1.5 eV. 

It is apparent from this figure that the decay rate drops 
dramatically as a result of the development of supercon- 
ducting order. From the calculation for model A, we see 
that there is a substantial drop resulting from the elim- 
ination of phase space. From model B, we see that the 
coherence factors reduce the decay rate further by a large 
amount. 

The absence of an antinodal quasiparticle peak in the 
pseudogap state and its subsequent emergence in the su- 
perconducting state has been interpreted here as the in- 
crease of its width as T c is approached. However, when 
the width becomes comparable to the quasiparticle en- 
ergy, i.e. as the curve reaches the dashed line in Figs. |31 
the quasiparticle concept breaks down. Once this occurs, 
our perturbative calculation can no longer be trusted, 
and it is not meaningful to assign a width or weight to 
the quasiparticle. What is significant is that it is possible 
to have a reasonably well-defined quasiparticle in the su- 
perconducting state as a result of phase space restrictions 



and coherence factors, as we have found. 

The broken symmetry state may or may not have well- 
defined quasiparticles at the single-particle gap edge. It 
depends on many non- universal details: the locus in mo- 
mentum space, the doping level, the interaction strength, 
etc. Thus, the absence of a well defined antinodal quasi- 
particle does not preclude a broken symmetry state. 
However, it may have important effects on non-universal 
aspects of the state, such as the temperature dependence 
of the order parameter which may, as a result, be strongly 
non-mean-ficld-likc. Also, our calculation leaves out fluc- 
tuation effects, which must be considered in the future. 

We end with three concluding remarks: (1) Although 
hole pockets cannot be observed in ARPES, other exper- 
imental probes can be used to look for their signature, 
for example, infrared Hall angle measurement in the un- 
derdoped regime^ (2) Because the interlayer tunneling 
matrix element is so strongly peaked at (tt, 0)^ we expect 
the c-axis optical conductivity to show a strong tempera- 
ture dependence at T c given our lifetime calculation. In- 
deed, this is consistent with the known measurements^ 
(3) We have not yet studied in detail the doping depen- 
dence. Nonetheless, it is possible to make a qualitative 
observation. There are two competing effects. As the 
system is moderately underdoped, the DDW order pa- 
rameter must increase. Thus, the quasiparticle in the 
(tt, 0) regime will have higher energy, increasing its scat- 
tering rate. On the other hand, the Fermi arcs will also 
shrink and the phase space for particle-hole scattering 
will decrease. We wish to return to these issues in the 
near future. 
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